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We review the recent developments in electronic Raman scattering measurements of charge ne¬ 
matic fluctuations in iron-based superconductors. A simple theoretical framework of a d- wave 
Pomeranchuk transition is proposed in order to capture the salient features of the spectra. We 
discuss the available Raman data in the normal state of 122 iron-based systems, particularly Co 
doped BaFe 2 As 2 , and we show that the low energy quasi-elastic peak, the extracted nematic suscep¬ 
tibility and the scattering rates are consistent with an electronic driven structural phase transition. 

In the superconducting state with a full gap the quasi-elastic peak transforms into a finite frequency 
nematic resonance, evidences for which are particularly strong in the electron doped systems. A 
crucial feature of the analysis is the fact that the electronic Raman signal is unaffected by the 
acoustic phonons. This makes Raman spectroscopy a unique probe of electronic nematicity. 
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I. INTRODUCTION. 


The study of correlated phases of matter obtained a big boost with the discovery of the iron-based superconductors 
(Fe SC) in 2008 CQ. These systems are interesting not just because they exhibit superconductivity at temperatures 
as high as 55 K, but also because they are a rich playground where the lattice and the electronic charge, spin and 
orbital degrees of freedom all play important roles. This complex interplay invariably leads to competition between 
various interesting phases that can be stabilized by varying temperature, by doping carriers and by applying pressure. 
Understanding this rich physics is a considerable challenge, and consequently the topic continues to be an active area 
of research in condensed matter systems HE]. 

Unlike the cuprate high temperature superconductors where the parent compounds are Mott insulators, the Fe 
SC systems are multi-band and multi-orbital metals at all doping. The undoped and the lightly doped compounds 
undergo a structural transition from a high temperature (T > Ts) tetragonal unit cell to a low temperature (T < Ts) 
orthorhombic phase (in the case of Fei+yTei-^Se^ the low-T phase is monoclinic for small x) which is followed in 
close proximity in temperature by a magnetic transition at T/v below which the system is antiferromagnetic. These 
transitions are suppressed with electron or hole doping, and beyond a certain doping value the system becomes su¬ 
perconducting with unusually high transition temperature Tc . Initial investigations of the superconductivity and its 
origin have focused mainly on the interplay between a stripe-like antiferromagnetism (where the magnetic ordering 
wave-vector is ( 7 r, 0) or (0, 7 r) in the IFe/cell notation) and superconductivity. While the issue is not entirely set¬ 
tled, it is popularly believed that the fluctuations associated with the stripe antiferromagnetism give rise to a s± 
superconducting pairing symmetry 0H6] in most, but possibly not all [7], Fe SC families. 

Ever since the reports of strongly anisotropic in-plane transport in the 122 systems [ 8 ] in the orthorhombic phase, 
a lot of attention has been given to study the property of nematicity in these materials. A nematic phase of matter 
is one which breaks rotational symmetry spontaneously, while preserving translational symmetry. Such phases have 
been studied extensively since the seventies in classical soft matter systems [Q], but relatively less is known about 
their quantum counterparts in electronic systems. However, their existence has been widely postulated for strongly 
correlated materials such as quantum Hall systems, cuprates, bilayer ruthenates mu, and most recently in the Fe SC. 
In the presence of a crystalline lattice, a nematic phase breaks discrete rotational symmetry, and as a consequence the 
associated order parameter is an Ising variable. In the context of the Fe SC, this order parameter is nonzero in the 
orthorhombic phase where the C 4 symmetry of the Fe unit cell is broken at the structural transition T 5 . Note that, 
in certain systems the structural and the magnetic transition are simultaneous ( Ts=Tn ), and, since the magnetic 
order by itself breaks C 4 symmetry, it is not clear whether the nematicity is a mere by-product of the magnetic 
order. Consequently, the issue of nematicity is more sharply posed for those systems where the structural transition 
precedes the magnetic one (Ts > T/v) leaving a finite temperature interval where C 4 symmetry is broken while the 
system remains paramagnetic P3HH]. The extreme example of this trend appears to be FeSe where only a structural 
transition [15] E>] is detected and the system remains paramagnetic down to its SC phase 021 CEE], hinting that nematic 
degrees of freedom may not be necessarily linked to magnetic ones. 

The microscopic origin of the nematic order is currently not known with certainty. One scenario is that the 
structural transition is, in fact, an instability driven by the anharmonic lattice potential, in which case the primary 
order parameter is the lattice orthorhombicity, and the electronic degrees of freedom are secondary order parameters 
that passively follow the symmetry breaking induced by the lattice strain. A second scenario is that the C 4 symmetry 
breaking is driven by electronic interactions, in which case the primary order parameter is electronic in origin. Within 
this picture one possibility is the spin-nematic transition whereby the spins of the two Fe sublattices phase-lock, 
which breaks C 4 symmetry, without developing a spontaneous magnetization, i.e., without breaking time reversal 
symmetry ESH2Z|. A second possibility is ferro-orbital ordering [28h32] . where either the occupations or the hopping 
matrix elements (or both) of the d xz and the d yz orbitals of Fe become inequivalent. Besides these two scenarios, other 
possibilities include a d-wave Pomeranchuk instability [33] in which the Fermi surfaces undergo symmetry breaking 
distortions due to interaction effects, as well as a valley density wave [34] . 

On the experimental side, initial studies have focused on the strong anisotropy of the electronic properties in 
the orthorhombic C 4 symmetry broken phase. Transport [ 8 l 1351 l39l . optical conductivity ]4QH44] . ARPES maim 
and neutron scattering j47h49] (reviewed in a separate contribution to this issue [50] ) performed on mechanically 
detwinned crystals all reported considerable electronic anisotropies. While it has been argued that the magnitudes of 
the measured anisotropies are too large to be due to the lattice orthorhombicity (which is 0.4 percent at most), such 
arguments can be at best quantitative, and therefore they do not convincingly rule out the lattice-driven scenario. 
Furthermore, even within the electronic-driven scenario, the above experiments cannot uniquely identify whether the 
primary order parameter is composed of electronic charge, spin or orbital degrees of freedom [51]. 

One difficulty in interpreting the above experiments is that, in the symmetry-broken phase, all the above order 
parameters are non-zero, and consequently it is difficult to identify the one which is most relevant. From this point 
of view it is desirable to set experiments that measure relevant susceptibilities in the C 4 symmetric phase, and search 
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for signatures of nematic fluctuations that soften upon approaching T#. One obvious possibility is the measurement 
of the orthorhombic elastic constant which measures the force constant associated with orthorhombic strain of the 
lattice [23l [52l [54] . Such studies are reviewed in a separate contribution to this issue [55]. However, note that an elastic 
constant measurement cannot, by itself, clearly distinguish between a lattice-driven from an electronic-driven scenario 
of nematicity. This is because an elastic constant measurement is a thermodynamic probe, and therefore, once the 
interaction between electrons with the acoustic phonons is taken into account, in both scenarios one would conclude 
softening of the relevant elastic constant. Next, in the spin-nematic scenario the fluctuations of the order parameter 
(which is a two-spin operator) involve a four-spin susceptibility. This poses a technical difficulty because, while NMR 
and neutron scattering can give crucial information on the evolution of the spin fluctuation spectrum HU EMU and 
the spin susceptibility at the antiferromagnetic wave-vector, they cannot easily access the four-spin susceptibility. This 
means that, at present there is no direct probe to test the fluctuations associated with spin-nematic order parameter. 
An alternative way to probe the nematic susceptibility was pioneered by Chu et al. [621464] who were able to extract 
the nematic component of the elasto-resistivity tensor in the tetragonal phase. This was achieved by performing strain 
dependent measurements of the transport anisotropy. Because the method allows a direct extraction of the lattice-free 
electronic nematic susceptibility, the observed divergence is a strong evidence for an electronic-driven structural phase 
transition. A drawback however, is the difficulty to associate the extracted nematic susceptibility to a microscopic 
nematic order parameter. 

On the other hand, in this review we argue that electronic Raman scattering allows a direct access to the fluctuations 
of the charge nematic order parameter, or equivalently the d-wave Pomeranchuk order parameter. This ability of 
Raman measurements has been somewhat overlooked in the past, although earlier Raman experiments in underdoped 
cuprates could possibly be interpreted along these lines [65]. We show in this review that Raman experiments in Fe 
SC give compelling evidences for the presence of nearly-critical charge nematic fluctuations in the tetragonal phase. 
These experiments also allow a direct extraction of the associated nematic susceptibility, which contains information 
about the incipient phase instability involving the purely electronic degrees of freedom. This is because the electronic 
Raman response is a spectroscopic probe, and it is “opaque” to the acoustic phonons in the system. In fact, this 
property of the electronic Raman response function provides a qualitative method to distinguish between lattice-driven 
versus electronically-driven scenarios of nematicity. It can be shown that in the former case the extracted nematic 
susceptibility from the Raman data should not show any signature of softening with lowering temperature. The fact 
that in Fe SC one does see softening proves conclusively that the nematicity is electron-driven rather than lattice-driven. 
This establishes Raman scattering as a key probe of translation symmetry preserving Fermi surface distortions which 
can be used to investigate other correlated electron systems such as the bilayer ruthenates and the cuprates where 
this kind of instability has been proposed but not yet confirmed unequivocally from an experimental point of view. 


This review is divided into two main parts, one theoretical and the other experimental. We start with the theory 
part which aims at giving a simple framework to understand how Raman scattering can be used to probe the charge 
nematic susceptibility and its associated dynamical fluctuations in an electron system. The approach is quite general, 
and the main features are expected to hold for the case of the Fe SC. In |II A| we consider a generic one band model 
with a charge nematic or Pomeranchuk instability where the Fermi surface breaks the C 4 symmetry. In |IIB 


provide analytical expressions for its critical fluctuations within the random-phase approximation (RPA). In |III A 


we 

we 


show that electronic Raman scattering directly couples to the charge nematic fluctuations, provided the appropriate 
symmetry channel is probed. In the case of the Fe SC this is the B\ g channel (which transforms as ( x 2 — y 2 ) 


the IFe/cell notation). We show in IIIB that the presence of these nematic fluctuations leads to the emergence of a 
quasi-elastic peak whose linewidth tends to vanish at the incipient pure electronic (ie, one without lattice coupling) 
nematic phase transition. We further show in |III Cj that Raman scattering measurements can access the charge nematic 
susceptibility, but only in the dynamical limit as opposed to the static limit which is relevant for the definition of 
the thermodynamic phase transition involving the C 4 symmetry breaking. While these two limits are the same for 
a pure electronic system, the coupling to the lattice induces a key difference between the two limits making Raman 
scattering measurement essentially blind to the linear coupling between the electron-nematic variable and the lattice 
orthorhombic strain. We conclude the theory part by discussing in |IIID| some of the additional subtleties associated 
with the multi-band nature of Fe SC. In particular we point out the existence of different flavors of charge nematicity 
when the orbital quantum number is taken into account. However, Raman measurements cannot distinguish between 
the various charge nematic order parameters that are possible in a multi-orbital environment. 

In the experimental part, after reviewing briefly the details of the Raman experiments in [TV[ we discuss the 
observation of charge nematic fluctuations in the tetragonal phase of electron-doped Ba(Fei_ a ,Co ;r As) 2 . In V A 


we focus on the behavior of the extracted charge nematic susceptibility as a function of Co electron doping and 
temperature. We then compare in |VB| Raman results with two other complementary probes of nematic fluctuations 
in the tetragonal state, elastoresistiviy and elastic constant measurements, and conclude that all three measurement 


are consistent with an electronic driven structural phase transition. We briefly comment on the role of disorder in V C 


by comparing Bal 22 and Sri 22 systems. In VIA we discuss the finite frequency spectra of the nematic fluctuations 
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and show that it is consistent with expectations from a simple mean-field approach of the nematic phase transition 
(VIB). We then show in VIC that Raman and shear modulus measurements can be consistently reproduced using 
a simple Landau-type mean-field picture with a linear coupling between charge and lattice nematicity. We conclude 
this review by addressing the fate of nematic fluctuations in the superconducting state (VII) and show that they can 
give arise to a novel collective mode, a nematic resonance, near the nematic quantum critical point. Experiments on 
electron doped Fe SC (Co-Bal 22 and Co-Nalll) appears to support its existence. 


II. THEORY: CHARGE NEMATIC INSTABILITY 

In this section we provide a microscopic description of charge nematic instability triggered by a phenomenologically- 
introduced electron-electron interaction. 


A. Model 

We consider a system of interacting electrons on a square lattice, simultaneously scattered by point-like impurity 
described by the Hamiltonian 


H = H 0 +H/ + V. (1) 

In the above 

^0 = e k c k ja c k,(7 (2) 

k,cr 

is the bare Hamiltonian of a band of electrons with dispersion ek, having lattice momentum k and spin a as quantum 
numbers, and described by usual creation/annihilation operators (c£ a , Ck,<r). Note that, while the Fe-based supercon¬ 
ductors (Fe SC) are multi-band systems, here, for the sake of simplicity, we restrict ourselves to a one-band model. 
Our main goal is to discuss certain qualitative physics, rather than quantitative ones, involving Raman spectroscopy 
near a charge nematic transition. We expect the main conclusions of this analysis to remain unchanged for a multi¬ 
band environment. Of course, in a multi-band system there are inter-band transitions that contribute to the Raman 
response, which are absent in a one-band model. However, typically such inter-band transitions are not related to 
criticality involving charge nematic transition. The latter is essentially a Fermi surface instability, and consequently, 
its critical fluctuations involve only intraband excitations. In Sec. EH we comment about the additional subtleties 
associated with nematic instabilities in a multi-orbital and multi-band environment. 

The interaction between the electrons is described by 

Ul = ( 3 ) 

q 


The operator 


On( q) 


VM 


£/*, q^k+q ,cr^k,cr i 


k ,cr 


( 4 ) 


where /k, q = (/ik + hk+ q )/2 with the B\ g (or equivalently x 2 — y 2 ) form factor = cos k x — cos k yi is the Fourier 
transform of the charge nematic operator (see Fig. [lj 



+ 


cJc i+£ + h.c.) 


-x-*y 


( 5 ) 


is the position of the lattice site i, and A f is the total number of sites. In the high temperature (^-symmetric phase 
(O n (r^)) = 0, while in the symmetry broken phase (O n (r^)) = 5o ? Vi, implying a charge nematic phase that preserves 
translation symmetry but not the discrete 7r/2 rotational symmetry, such that the hopping matrix elements along x 
and y are inequivalent. This is the lattice version of a Pomeranchuk transition P33E5I where a spherical (or circular 
in two dimensions) Fermi surface becomes ellipsoidal (or elliptical), and indeed in the following we do not distinguish 
between the names “charge nematic” and “Pomeranchuk”. We take the constant go > 0, implying that C 4 symmetry¬ 
breaking distortions of the Fermi surface are energetically favored by the interaction term. This is a pre-requisite for 
a charge nematic instability, at least in weak coupling theories involving random phase approximation. 
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FIG. 1. Left: B\ g symmetry form factor hk of the Pomeranchuk or charge nematic operator in k-space. Right: associated 
Fermi surface deformation in a one-band model. The original spherical Fermi surface is distorted into an elliptical one in the 
charge nematic phase. 


Evidently, there is non-trivial physics involved in transforming the usual short-range repulsion between the electrons 
into an interaction given by Eq. <§, and in the Fe SC this can be a consequence of spin, charge, or orbital fluctuations 
[SU (The nature of magnetic interactions and their interplay with orbital degrees of freedom in Fe SC is reviewed in a 
separate contribution to this issue El)- However, an inquiry into the origin of this interaction is related to the question 
as to what triggers the nematic transition. To the best of our knowledge, this issue is not entirely settled for all the 
Fe-based systems, and is beyond the scope of the current review. Instead, here we adopt a more phenomenological 
point of view whereby, assuming the existence of such a transition, we write down a minimal effective theory that 
describes the transition and the critical phenomenon associated with it. It is in this low-energy effective theory sense 
that one should understand the above interaction term. 

The last term 


^ = ^oEL q /M, (6) 

k,q,cr 

describes scattering of electrons with isotropic point-like impurity potentials. The role of this term is to provide 
a finite lifetime to the electrons. In the following the impurity term does not affect the description of the charge 
nematic instability. However, the inclusion of finite lifetime is crucial for a meaningful discussion of the Raman 
response function. In a perfect metal, where quasi-particles are infinitely long-lived, the Raman response, involving 
intra-band particle-hole excitations and zero momentum transfer with the photons, vanishes. This is due to the fact 
that the constraints from energy and momentum conservation cannot be satisfied simultaneously. Adding impurity 
scattering is an effective way to bypass the latter, and to take into account the intra-band contribution that is always, 
in practice, present. 


B. Instability and critical fluctuations 

We first discuss the effect of the impurity scattering. We assume that the impurity potential is weak enough such 
that it can be treated in Born approximation. This provides a finite lifetime to the electrons given by 

t -1 = 27mV 0 2 Po, (7) 

where rii is the impurity concentration, and po is the density of states at Fermi energy. 

Next, we treat the interaction in random phase approximation. The description of the charge nematic instability 
is facilitated by introducing the Hubbard-Stratanovich field </> n (q) to decouple the interaction which can be rewritten 
as 


Hi = y F + 0n(-q)O„(q)]. 

q 


( 8 ) 
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The second term above describes interaction between electrons and the Hubbard-Stratanovich field 0 n (q), which is 
shown graphically in Fig. [ 2 ] (a). With this rearrangement the theory is formally quadratic in the fermionic variables, 
which can be integrated out. This leads to the action in terms of the critical variable 

S[pn] = Wn)\ 2 , (9) 

q ,iv n 

where % n (q, iv n ) is the nematic susceptibility given by 

Xn^q^n) =fl , o[l-fl'on„(q,iz/„)], ( 10 ) 


with the nematic polarization 


n„(q,*z/„) = - 


2 

P 


^ [ /k,q^k(^n)^k+q(^n 
o; n ,k 


(ii) 


In the above /3 = 1 /(fc^T), where kp is the Boltzmann constant, and the factor 2 is due to summation over the spin 
index. The electron Green’s function is given by 

G k (iw„) -1 = iu n ~ e k + */(2r)Sgn(w n ). (12) 

Note that, since the impurity potential is isotropic, for symmetry reasons the vertex correction due to impurity 
scattering affects neither the q —>> 0 limit of Xn(q 5 0), which controls the charge nematic transition temperature, 
nor the susceptibility x n (0,iz/ n ), which will be relevant when we discuss the Raman response in the next section. 
Consequently, in the current model the role of the impurity is restricted to providing finite lifetime to the electronic 
excitations. 

It is convenient to divide the frequency and momentum dependencies of the nematic polarization into two parts, 
namely the contributions from the high-energy and the low-energy electrons, such that 

nn(*T Wn) n n (0, 0) T n n (q, il^n )high T n n (q, il^n )low (13) 

Since we do not expect any singular contribution from the high-energy electrons, the resulting frequency and momen¬ 
tum dependence is analytic. To lowest order we expect 

n„(q,iz/ n )hi g h = ci(q/k F ) 2 + c 2 (v n /E F ) 2 . (14) 

In the above II n (0,0) and the coefficients depend on the details of the band structure. The charge nematic 
transition takes place when the Stoner criterion 

ro = l-^ 0 n n (0,0) = 0 (15) 

is satisfied. In general ro, which is related to the nematic correlation length £ n by ro = (a/in) 2 with a being the unit 
cell length, is temperature dependent and it decreases as the nematic correlation length increases upon approaching the 
instability with lowering temperature. In electron-doped Bal22 it is now known from elastic constant measurements 
(discussed in a separate contribution to this issue [55]) and also from Raman scattering measurements (described later 
in this review) that ro is linear in temperature over a wide range, i.e. 

r 0 (T)=fo(T-T 0 ), (16) 


where Tq is the charge nematic transition temperature. 

The low-energy contribution II n (q, iv n )\ ow determines the dynamical properties of the charge nematic fluctuations, 
and its evaluation is quite standard. The ek-integral can be performed by linearizing the electronic dispersion. This 
gives to leading order 


(Qj ^n)low i'L'nPO 


p 27 T 

Jo 


d 0 k 


hi 


2 i r iS Un — r j pgcos(6 > k — 0 q ) 


In the above S Uri = u n + Sgn(z/ n )/r, vp is Fermi velocity, and 6\< is the angle of k measured from one of the two 
equivalent major axes of the unit cell. The angular integral can be performed analytically if we approximate the B\ g 
form factor by « — cos(2#k)- From the angular integral we get 


n„(q,iz/ n )low = [cos 2 (2#q) 1(7(flq,!/„) 

V pQ 

+ sin 2 (20 q )/,s(a qiJ/ J] , 


(17) 




where a qj „ n = [v n + Sgn (u n )/T]/(v F q), 


Ic(a ) 



d9 cos 2 29 
27r ia — cos 0 


i(l H- 2a 2 ) 


(1 + 2a 2 ) 

_ \/l + a 2 


Sgn(a) — 2a 
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(18) 


and 



dO sin 2 2 6 
27r ia — cos # 


= 2ai 


1 2a 2 


2 |a| \/l + a 2 


(19) 


Thus, the overall q-dependence of the nematic polarization is anisotropic, which is eventually a consequence of the 
form factor associated with the nematic variable defined in Eq. Q [68] • In the above the ratio a q? ^ n is large for 
temporal fluctuations and it is small for spatial fluctuations. 

Quasi-static limit, a q ^ n 1: This limit is relevant for studying the thermodynamic signatures of the charge 
nematic instability. Using the properties 


Ic(cl —» 0) = iSgn(a) — 2ia, 
Is(a — y 0) = 2ia, 


we get in this limit 


Xn (q,*Z / n)=50 


- 4 c 3 


r 0 + Ci JTT + C -2 + 2C3 m cos 2 2# q 




{y F q) [ 


r(l + i\ v n\r) 1 ) cos4$, 


( 20 ) 


with c 3 = gop 0 . 

Quasi-dynamical limit, a q ^ n 1: This limit is relevant for studying the signatures of the critical mode </> n (q) in 
Raman spectroscopy [69 . This is because in typical Raman scattering experiments the momentum q given to the 
electrons by the visible photons is negligible compared to the Fermi momentum. Using 


/c(a oo) = Is(a oo) = i/(2a), 


we get in this limit 


Xn^Q ,i”n) =90 


q 2 , ^ , K| 

ro+Cl g: + ‘ :2 £|: +C3 wTT7Tj' 


( 21 ) 


In this limit the momentum anisotropy is absent from the leading order terms. 

Analyticity of Xn{ 0,0): From Eqs. (20) and (21) it is important to note that y n (0,0) obtained from the static and 
the dynamic limits are the same, i.e., for an infinitesimal g > 0 


lim lim y n (q, uj + ip) = lim lim y n (q, uj + ig) = -. 

g —^0 cj —)-0 u —^-0 q —±0 gtTo 


( 22 ) 


This implies that y n (q, a; + iv) * s analytic at zero momentum and frequency. 

As a quick aside, the above behavior is to be contrasted with the susceptibility of a conserved quantity such as 
the charge susceptibility, which as we will show below is not what is measured in the Raman response. The charge 
susceptibility is defined by 

Xc(q,Wn)=/ (T r p_ q (r)p q (0))e* ly ” T , (23) 

J 0 

where T r is the imaginary time ordering operator, and 




k+q ,a^k,o 


k,cr 


(24) 
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FIG. 2. (a) Interaction between electrons (solid lines) and the charge nematic boson (wavy line) (j) n . The interaction vertex 

has nontrivial momentum dependence given by /k, q = (hk + hk+ q )/2 with the B\ g form factor hk = cos k x — cos k y . (b) Graphs 
for the Raman response function. The first term is the quasiparticle contribution which is non-critical. The second term is 
the contribution of the nematic boson which has information about criticality associated with the nematic instability. This 
contribution is non-zero only in the B\ g Raman channel. 


is the Fourier component of the charge density operator. In this case it is well known that the uniform charge 
susceptibility from the static limit is finite with 

lim Xc(q,0) = Xc 7 ^ 0, (25) 

q—m 

where %c is the charge compressibility of the electronic system. On the other hand the dynamical limit vanishes with 

Xc(0,w + ir?) = 0. (26) 

This is a consequence of particle number conservation, i.e., [H,p q =o] = 0, which in turn follows from the global 
U{ 1) symmetry of the Hamiltonian. In other words, this symmetry enforces a non-analyticity at zero frequency and 
momentum, and 

lim lim Xc(q 7 ^ + irj) ^ lim lim Xc(q?^ + ir}). (27) 

<?—>-o to —m id —m 0 


By contrast, the uniform nematic operator is not a conserved quantity, since ['H/,O n (q = 0)] ^ 0, and therefore 
there is no physical reason to expect a similar non-analyticity in Xn (q, w + irj) at zero frequency and momentum in 
purely electronic models , i.e., those where the coupling of the electrons to the lattice strains is ignored. As we discuss in 
the next section, the analyticity implied by Eq. ( [22] ) is important for interpreting the signature of the charge nematic 
instability in Raman response. 

Note that, in order to establish the analyticity of Xn(0,0) it is crucial to consider electrons with finite lifetime. It 
is easy to check that Eq. (22) does not hold for an ideal metal for which r _1 — 0. However, such a non-analyticity, 
which is not associated with any symmetry, and whose origin can be traced to the fact that for an ideal metal the 
phase space for particle-hole excitations is sharply defined, is rather an artefact. In practice, the electrons have a 
finite lifetime, and this ensures that the phase space for particle-hole excitations is no longer sharply defined. 


III. THEORY: ELECTRONIC RAMAN RESPONSE NEAR CHARGE NEMATIC INSTABILITY 

In this section we discuss the characteristic signatures of a charge nematic instability in electronic Raman spec¬ 
troscopy. We also discuss how electron-lattice coupling affects such an instability. We argue that the Raman response 
is “blind” to this coupling, and therefore it is an ideal tool for studying the bare electronic nematic correlations. 


A. Big response, static and dynamic limits 


The theory underlying the electronic Raman spectroscopy for correlated systems has been reviewed elsewhere 
[221 E]. Here, we simply remind few salient points associated with this experimental technique. Accordingly, we 
define the stress tensor by 


T a/3 (q) = E 

k,cr 


<9 2 e k f 

^fc a ^fc/ k +‘l/ 2 ’ ffCk - q / 2 ’ <J, 


( 28 ) 
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and the associated correlator as 


Afl n r 


Xa/3(q,ifin)= [ dT(r T rL(q,T)r a/3 (q,0))f 

Jo 

Following Kubo formalism, the analytic continuation of the above gives the response function 


(29) 


Xap(<l,fy 


lim Xap(qJ^n) 

rOO 

i / dte lM { T^(q,i) , T a0 ( q,0) 
Jo 1 J 


(30) 


The imaginary part of the stress tensor response function at its dynamical limit, i.e., X^si *4 = 0, H) = Imy a ^(q = 0,12), 
is accessible by means of Raman scattering experiment. This is because the associated scattering cross-section is 
proportional to the correlation function S a p(Q) = (T^p(Q)T a p(Q)), which in turn satisfies the fluctuation-dissipation 
theorem 


s a p(ft) = 1 [1 + n B (fi)] Xa/3(q = fi )> 


(31) 


where tib(I 2) is the Bose function. 

In the Fe SC, using the notations of a unit cell with IFe/cell, the quantity of interest is the B\ g stress tensor given 
by 


^B lg (q) — T xx (s\) Tyyio)- (32) 

Within the effective mass approximation [72] . valid for non-resonant Raman scattering, we define the associated 
Raman vertex as 


7 s l 9 (k) = (pl x - djQ e k = 2*i h k , (33) 

where the dispersion of Eq. Q is given by ek = —2ti(cos k x + cos k y ) H-, with t\ being the nearest neighbor hopping 

parameter. In principle, k) can include higher harmonics of the same B\ g symmetry, such as that coming from 
the third nearest neighbor hopping, if present in J-Lq. In the following we ignore such terms since they do not affect 
the results qualitatively. 

The computation of the B\ g correlator within random phase approximation (RPA) is quite straightforward. We 
get 


Xs l9 (q,*^n) = n 77 (q,*n ra ) + 5 0 n 7 ^(q,*O n )x„(q,iO n ), 


(34) 


where II 77 is defined like II, 


in Eq. (11) with the form factors q replaced by 7 # lfl (k), and ITy^ is defined similarly 
lfl (k)/ik- Tlie graphical representation of these two terms is shown in Fig. [21(b). Note that, 


with the form factors 7 b 

the second term, which contains the critical contribution, is non-zero only in the Raman B\ g channel! In the A\ g and 
the B^ g channels this term is zero by symmetry. Indeed, in experiments the signature of the charge nematic instability 


is observed only in the B\ g Raman channel. Since, II 77 = At^U. n 
momentum 


and n 7/l = 2 tiII n , we get, at low frequency and 


XBig (Q: i^n) ~ 4t 1 % n (q, i!2 n ). (35) 

The above proportionality implies that the properties of y n (q, zI 2 n ), discussed in Sec. [TIJ are also relevant for 
X#i fl (q, ^n) (see also [3, El])- In particular, in purely electronic models the B\ g response function is analytic 
at zero frequency and momentum. 

In electronic Raman spectroscopy in the B\ g geometry we measure a quantity proportional to Xb 19 (q = 0, II) = 
Imy^ lg (q = 0,11). This quantity can be used to deduce the frequency-integrated spectral weight of the associated 
Raman conductivity Xb 19 ( q = 0 ,12)/I2 which, by Kramers-Kronig relation, gives the real part of the uniform response 
function. Thus, 


dynamic _ ^ 

= n 


L 


00 dSl .. , 

-fl-XB lB (q = o,n) 


= i im „XBi 9 (q = o^)' 
0^0 y 


( 36 ) 
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It is important to distinguish the two conceptually distinct quantities X^ amic and 


static _ 

x Blg - 


l™ XB lg (q, = 0). 

q^O y 


(37) 


Note that, the thermodynamic instability at a charge nematic transition is associated with a divergence in the latter 
quantity. However, since XB lg (q, O) is analytic at zero frequency and momentum in purely electronic systems, we get 

xt n ; mic = xtl iC ’ ( 38 ) 

and, therefore, the divergence is “visible” from the dynamical limit. In other words, one can obtain information about 
the divergence of a susceptibility associated with the charge nematic instability from Raman spectroscopy. 


B. Signature of instability: quasi-elastic peak 


In the following we study the low frequency properties of the Raman response near the transition. Using Eqs. (21) 
and (35), for U < 1/r <C Ep we get 


xsi„(q = o,n) = A 0 


ro + c 3 - 


o 


1 -1 


U + i/r 

with Aq = At\/go and C 3 = gopo • This implies a Raman response with 


x' Blg ( q = o,fi) = 


-1 


n 


a 0 c 3 t 
(r 0 + c 3 ) 2 V ^ 2 + r 2 


where 


r = 


ro 


(ro + c 3 ) 7 


(39) 


(40) 


(41) 


Thus, the characteristic signature of the transition in the low frequency Raman conductivity is a quasi-elastic peak 
with Lorentzian lineshape that sharpens as the system approaches the transition, since the width T —>• 0. 

It is useful to note that the deduced from the experimental data [75] (described in |v| has the form 


dynamic _ ^0 


XB! 


r o 


(42) 


where B is a non-singular part that is often temperature independent. This apparent violation of the Kramers-Kronig 
relation is partly due to the fact that, in practice, the upper cutoff of the frequency integral of Eq. (36) is finite and 
is set to a value beyond which the measured Raman spectra is temperature independent. A second reason for the 
H-term is that part of the electronic Raman signal observed is symmetry independent and therefore unrelated to 
criticality. 

Finally, we note that, besides Xj 3 ^ amiC an( ^ width T, a third quantity which is experimentally accessible is the 
slope of the Raman response at zero frequency 


Si m 


dnXB lg ( q = °,^) 


A 0 c z t 


J 0=0 


(43) 


where the equality is the result of the random phase approximation. Using T ~ ro/(c 3 r) close enough to the phase 
transition, we get a scaling relation between the three experimentally accessible quantities 


dynamic 
XB\ g 


r s t , 


(44) 


provided we ignore the non-singular H-term in Eq. (42). This provides an independent check for the validity of the 
random phase approximation theory. 
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C. Effects of coupling to the lattice 

A crucial ingredient in the above discussion is the analyticity of the nematic susceptibility y n (q, iH n ), and therefore 
that of the correlator Xs lg (q, i£l n ) involving the B\ g stress tensor, at zero momentum and frequency. It is this 
property that allows us to conclude that the thermodynamic divergence at the charge nematic phase transition, which 
shows up in the susceptibility taken to its static limit, is also “visible” from the dynamic limit which is accessible 
via Raman spectroscopy. As we discussed earlier, in order to demonstrate the analyticity, it is important to consider 
electrons having a finite, frequency independent, lifetime which, in simplest models, can be typically attributed to 
impurity scattering. 

The above conclusion, however, holds only if we ignore the symmetry-allowed coupling of the electronic charge 
nematic operator O n ( q) of Eq. © to the orthorhombic strain of the underlying lattice. In practice, however, such 
electron-lattice coupling is always present in crystalline solids. In the Fe SC, the fact that the C 4 symmetry breaking 
(ignoring the positions of the As and the chalcogen atoms) is invariably accompanied by orthorhombic distortion, 
provides phenomenological proof of the presence of such a coupling. Consequently, it is worthwhile to examine how the 
above statements concerning the charge nematic transition and its Raman “visibility” are modified once the coupling 
to the elastic strain is taken into account. 

In the following we consider a two dimensional square lattice whose elastic free energy, to lowest order in the strains, 
is given by 


Fe — (e xx + €yy) + \ y + Ci 2 e xx e yy . (45) 

Here e^- = {diUj + djUi)/ 2 , with (i,j) = (x,y) are the strains, the vector u denotes displacement from equilibrium, 
and C 11 etc. denote elastic constants in Voigt notation. 

In terms of the above the orthorhombic strain is given by es( r) = e xx (r) — e yy {r), and the accompanying bare elastic 
constant is C$ = (Cn — C i 2 )/ 2 . We write the symmetry-allowed coupling between the orthorhombic strain and the 
electronic charge nematic operator O n ( q) as 


Hel —lattice — ^0 y]On(q)es(q), 

q 


(46) 


where es{ q) is the Fourier transform of es (r), and Ao is the coupling constant having the dimension of energy. 

Studying the detailed implications of the above coupling on an electronic nematic phase transition is beyond the 
scope of the current review, and will be presented elsewhere. Here we focus only on the following salient points that 
are relevant to the current discussion. 

(i) One effect of the coupling is to increase the temperature of the C 4 symmetry breaking nematic/orthorhombic 
transition from To, defined in Eq. (16), to 


T S = T 0 + X 2 0 /(C° s f 0 g 0 ). 


(47) 


Below Ts the C 4 symmetry breaking is manifested both in the electronic sector, where the dispersions becomes C 2 
symmetric, as well as in the lattice sector with orthorhombicity es 7 ^ 0. Note that the effective orthorhombic elastic 
constant can be expressed as (for a derivation see, e.g., Ref. [55]) 


C s = C° s - A§x n (0,0), 


(48) 


where Xn( q,^) is the bare electronic nematic susceptibility defined in Eq. (10). This implies that Cs vanishes at 


T = Ts, and consequently the transition can be detected in experiments that measure Cs either directly by ultrasound 
or indirectly by bending techniques. 

(ii) At Ts only the effective orthorhombic elastic constant Cs vanishes, while the remaining elastic constants stay 
finite. An important consequence of this is that the critical fluctuations are restricted to two high-symmetry lines 
q x = ±q y in the two-dimensional Brillouin zone [76, [77]. This can be understood from the following. Writing 
q\ = ( q x + q y )/V 2 , it can be shown that along the lines q x = q y and for the polarizations ni = ( 1 ,- 1 ) the acoustic 
phonon dispersion is given by Hi q == {Cs / pY^qi, where p is the atomic mass density. Similarly, along q x = —q y and 
for the polarization n 2 = (1,1) the dispersion is H 2q = (Cs/p) 1 ' 2 ^, where q 2 = {q x — q y )/y/2. It can be shown that 
these are the only two directions in the Brillouin zone for which the phonon velocity vanishes at Ts- This is because 
for all other q the acoustic phonons excite not just the critical strain es, but also the non-critical ones whose elastic 
constants remain finite at Ts. 
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(b) 



d-wave Pomeranchuk 


(c) 0 



xz/yz ferro-orbital order 


FIG. 3. (a) Typical Fermi surface of most Fe SC within a 3 orbital tight binding model (see e. g. Daghofer et al. [79]). The 

orbital content of each pocket is indicated. The B\ g form factor cos k x -cos k y is depicted in the background, (b) Fermi surface 
deformation associated to the d- wave Pomeranchuk order Oq (c) Fermi surface deformation associated with xz/yz ferro-orbital 
order O s Q with gk — 1 [28l |29j . 


(iii) Finally, and most importantly for the current discussion, as a result of the electron-lattice coupling the nematic 
susceptibility acquires a non-analytic correction. Denoting the dressed susceptibility as % n , for q = q a q a (summation 
not implied) with a = (1, 2), we get that 


Xn 1 (« laQcnittn) = Xn 1 (QaQa, *^n) ~ 


CWa + W 


(49) 


Here q a are the unit vectors along the two critical directions within the Brillouin zone. The non-analyticity of the 
second term is a consequence of translation symmetry. Since moving all the atomic positions by a fixed displacement 
does not change the overall energy of the system, the electron phonon coupling, which is the numerator of the second 
term above, vanishes in the uniform limit q a 0. As a consequence, the effect of the acoustic phonons is entirely 
absent in the dynamical limit and x n (q = 0, fi) = x n ((q = 0,0)) [78] . In other words, the Raman response, being 
opaque to the acoustic phonons, measures the bare electronic nematicity and X^ amic oc (T — To) -1 deduced from 
it tends to diverge at the purely electronic temperature scale Tq [69]. This is in contrast to the inverse of the 
orthorhombic elastic constant C^ 1 oc (T — T^) -1 . However, the divergence of Xb™™ 10 * s invariably cutoff at the 
actual C4 symmetry breaking transition Tg. Conversely, since the electronic Raman response function is unaffected 
by the acoustic phonons, any signature of nematicity seen in Raman is unambiguous proof that it is electronic in 
origin. 


D. Extension to multi-orbital systems 

Until now the qualitative physics was described in terms of a single band system in order to simplify the discussion. 
Nevertheless, it applies equally well to multi-orbital (and multi-band) systems like the Fe SC. In a multi-orbital 
environment the main additional novelty is that, due to the presence of the orbital quantum number, one can construct 
different flavors of charge nematic order parameters. The most direct extension of the single band case described in 
Eq. @ is a multi-orbital version of the d-wave Pomeranchuk instability where the orbitally resolved electron densities 
have all B\ g nematic form factors : 


O 


d _ 
o — 


1 

77 




(50) 


where n k is the electron density operator, m is the orbital index, and h k = cos k x — cos k y . While little discussed 
initially this instability has recently been put forward as a candidate order parameter for the orthorhombic phase of 
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FeSe [8QM82] . A second possibility is to define a orbital dependent order parameter of the form : 


O 


s 

o 


1 

VN 






(51) 


involving the xz and the yz orbitals. Here g k is a function with A\ g symmetry such as = 1 (equivalent to xz/yz 
ferro-orbital order [28H3Q1 183] ). or g k = cos k x +cos k y . Note that, the actual deformations of the Fermi pockets in the 
nematic phase, and the axis along which a pocket will elongate/contract, can vary depending on the choice of g^. In a 
two-orbital model relevant for Fe SC the deformations obtained with g^ = 1 are qualitatively similar to that obtained 
with a non-zero d-wave Pomeranchuk order parameter 0%. One reason for this is that the electron pockets composed 
mostly from the xz and the yz orbitals are centered around (7r, 0) and (0,7r) respectively, and the form factor dk is 
approximately a constant, with opposite signs for these pockets. Consequently, the projections of these two order 
parameters on the electron pockets are indistinguishable. The various Fermi surface deformations associated with 
different choices of the nematic order parameter are illustrated in Fig. [3j 

Overall, we notice that in a multi-orbital system the nematic order parameter can have non-trivial structure in the 
orbital space, and it can also be accompanied by momentum space structures that are different than what is possible 
in the single-band case described by Eq. From the point of view of Raman spectroscopy it is important to note 
that, irrespective of the details of their momentum and the orbital space structures, as long as the nematic order 
parameter transforms as a B\ g object, the critical fluctuations will be observable in the electronic Raman B\ g channel. 
This can be illustrated by calculating the form of the B\ g stress tensor or vertex for a specific tight binding model of 
the Fe SC. In a multi-orbital system one can generalize Eq. (33) and write the component of the non-resonant B\ g 
Raman vertex in the effective mass approximation |84] 


TBi" = (d] 




(52) 


where e 
(m, n=d. 


rnn 

k 


are the components of the tight binding dispersion in the orbital basis. In the minimal two-orbital 
, d yz ) model of the Fe SC of Raghu et al. [55], the B\ g vertex is diagonal in orbital space with 




(53) 


where t\ and 12 are the near-neighbor hopping parameters for a and tt type Fe orbital overlap respectively. We note 
that next nearest-neighbor hopping integral along the diagonals of the Fe square plane do not contribute to the B\ g 
Raman vertex. The k-space structure of the two diagonal Raman vertex matrix elements and the associated Fermi 
surface deformation is illustrated in Fig. [4] for the tight binding parameters of Raghu et al. [85. 

In this model the B\ g stress tensor is 

Tb 19 (q = 0) = (ti + t 2 )O d 0 + (h - t 2 )O s 0 , (54) 

with $k = cos k x + cos k y . The B\ g stress tensor has finite overlap with both the ferro-orbital and the d- wave 
Pomeranchuk order parameters, and therefore criticality in either of these two channels is manifested in the Raman 
Bi g response. Conversely, based on the Raman data, it is difficult to determine whether the nematic criticality is 
ferro-orbital or d- wave Pomeranchuk type. 


E. Summary of the theoretical discussion 

Here we summarize the main points of the theory developed in the previous and in the current sections. In Sec. [H] 
we developed a Drude-RPA theory to describe a charge nematic or d-wave Pomeranchuk phase transition, starting 
from a phenomenological electron-electron interaction. In particular, we showed that in purely electronic models 
the nematic susceptibility y n (q, fi), that describes static and dynamic fluctuations of the charge nematic operator 
O n ( q), is analytic at zero momentum and frequency, provided we take into account the effect of impurity induced 
elastic scattering of the electrons on their lifetimes. The latter ensures that the single particle scattering rate does 
not vanish at zero frequency, which is a crucial ingredient in order to establish the analyticity. We also noted that 
this analyticity is eventually tied to the fact that the uniform charge nematic operator O n (q = 0) is not a conserved 
quantity. Next, in Sec. |TTT] we discussed the characteristic signatures of the charge nematic instability in the electronic 
Raman response. Within the Drude-RPA theory the effect of the criticality is symmetry-selective in the sense that it 
is observed only in the Raman B\ g channel, and not in the other Raman channels. We showed that the B\ g Raman 
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FIG. 4. (a) k-space structure of the B\ g Raman vertex matrix elements in the orbital basis for the two orbital model of Raghu 
et al. m (see also m for a more realistic 5-band calculation of the Raman vertex). In this simplified model the main features 
of the Fermi surface topology of Fe SC can be reproduced with t 2 —-l.Sti. (b) Fermi surface deformation associated with the 
corresponding B\ g average stress tensor projected at the Fermi wave-vectors within the 2 orbital model. 


response (fi) is proportional to the imaginary part of the nematic susceptibility in its dynamical limit, i.e., to 
Imx n (q = 0, fi). Next, using the analyticity discussed above, as well as the Kramers-Kronig relation, we argued that 
the frequency-integrated Raman conductivity Xs lfl (^)/^ i s a measure of the nematic susceptibility at its static limit, 
i.e., liuiq-^o Xn(<F ^ = 0), which is the quantity whose divergence signals the second order charge nematic transition. 
Furthermore, we showed that when the single particle lifetime is dominated by elastic scattering, the B\ g Raman 
conductivity has a Lorentzian lineshape whose width narrows as a function of temperature as the system approaches 
the nematic instability. Finally, we pointed out that the electronic Raman response function naturally screens out the 
effect of the coupling of the electronic nematic variable with the orthorhombic strain of the lattice. Consequently, the 
Raman response is a measure of the bare electronic nematicity that is unaffected by the presence of the lattice. As 
such, it is an ideal tool for providing qualitative distinction between an electronically-driven nematic instability from 
a lattice-driven one. 


IV. RAMAN EXPERIMENTS 


Raman scattering is a photon-in photon-out process in which an incident photon with energy ujl and polarization 
cl is inelastically scattered by the medium into a photon of energy ujs and polarization es (Fig. [5|. For a Stokes 
process, uos < ojl, an excitation with energy uj = ujl — us is created in the solid, u) is usually referred to as the 
Raman shift and is traditionally given in units of cm -1 (8.066 cm -1 =l meV). A typical Raman spectrum for a metal 
consists of sharp peaks due to Raman allowed optical phonons superimposed on a continuum of electronic origin. 
The Raman experiments described here were carried out using a diode-pumped solid state laser emitting at 532 nm 
or a Ar-Kr mixed gas laser with several lines in the visible spectral range. The inelastically scattered photon were 
analyzed using a triple grating spectrometer equipped with a nitrogen cooled CCD camera. Special care was taken 
in order to determine the laser induced heating. It was first estimated by comparing the power and temperature 
dependencies of the phonon frequencies. This estimate was then cross-checked by monitoring the onset of Rayleigh 
scattering by orthorhombic structural domains across the structural transition temperature as a function of laser 
power. For Co-Bal22 crystals both methods yielded an estimated heating of 1 K db 0.2 per mW of incident power. 
In order to extract the imaginary part of the Raman response function, the raw spectra were corrected for the Bose 
factor using Eq. (31) and for the instrumental spectral response as well. 


In this review we will mostly discuss the spectra performed in the B\ g symmetry because they probe the nematic 
degrees of freedom relevant to the Fe SC. The x 2 — y 2 or B\ g symmetry can be selected by choosing crossed incoming 
and outgoing photon polarizations at 45 degrees with respect to the Fe-Fe bonds. Here the notation B\ g refers 
to the one Fe unit cell whose axes are along the Fe-Fe bonds. For comparison we will also show spectra in B^ g 
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photon out 
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B 2g / xy 



Bi g / x 2 -y 2 

• • 



FIG. 5. Sketch of the scattering geometry in a Raman experiment. The photon polarization configuration corresponding to 
the Big and B<i g symmetries are depicted with respect to the FeAs plane. Note that the 1 Fe unit cell notation is used. 


symmetry which can be selected by choosing incoming and outgoing photon polarizations along the Fe-Fe bonds. 
Switching between both symmetries was usually performed by rotating the crystal by 45 degrees while keeping both 
the polarizer and the analyzer fixed. Note that in terms of the full lattice unit cell (or 2 Fe unit cell), which has its 
axes at 45 degrees to the Fe-Fe bonds and is sometimes used in the literature, the B\ g (B 2 9 ) symmetry discussed here 
corresponds to the B^g (B\ g ) symmetry. 


V. EXPERIMENTS: CHARGE NEMATIC SUSCEPTIBILITY OF ELECTRON DOPED 122 SYSTEMS 

A. Electron doped Ba(Fei_ x Co x As )2 


We first discuss Raman data in the normal state of electron doped Ba(Fei_ a ,Co ;E As )2 (Co-Bal22) and the evolution 
of the extracted static nematic susceptibility as a function of temperature and Co doping [75] . Being one of the 
most studied Fe SC system, the salient features of the phase diagram as a function of Co doping are well-known [86]. 
There is a quasi-simultaneous magneto-structural transition in the parent compound (Ts ~ T/v) [87, §8] which splits 
under Co doping with Ts > T/v HD Q2J HMSt The superconducting dome starts at x ~ 0.02 and extends up to 
at least £=0.15. All the Co contents x mentioned in this review were determined using wavelength dispersing X-ray 
spectroscopy (WDS). 

The low energy (uj < 600 cm -1 or 75 meV) Raman responses of the parent compound BaFe 2 As 2 in the B\ g and 
B^g symmetries are shown as a function of temperature in Fig. pna) and (b). In the tetragonal phase the B\ g response 
shows a strong enhancement upon cooling towards Ts =138 K, before collapsing in the orthorhombic / spin density 
wave (SDW) state. By contrast the B^g response is essentially independent of temperature in the tetragonal phase 
and only shows a mild suppression below Ts. The observed symmetry dependence is in agreement with the presence 
of dynamical nematic fluctuations having x 2 -y 2 symmetry as discussed above. This interpretation is confirmed by 
the Co doping dependence of the spectra which shows a systematic enhancement of the B\ g response towards Ts and 
the disappearance of any temperature dependence in the strongly electron overdoped composition, far away from the 
orthorhombic instability (£=0.20, see Fig. [7]). 

From the behavior of the Raman response at finite energy one can define the associated nematic susceptibility 
where /i is the symmetry channel, using Kramers-Kronig relation (see Eq. 36). As emphasized above, the 
susceptibility is in this case obtained in the dynamical limit. The physical quantity governing the nematic susceptibility 
in this limit, the Raman conductivity x^/cj, is shown as a function of symmetry for Bal22 in Fig. [6j It is also shown 
as function of electron Co doping in Fig. [7] While the Raman conductivity is flat in B 2 g symmetry, in B\ g symmetry 
it is dominated by a peak centered at zero energy over a wide range of Co doping in the tetragonal phase. For 
£ <0.045 the peak amplitude grows upon approaching Ts and collapses quickly below. It is interesting to note that 
the enhancement of the peak amplitude is also seen down to the superconducting phase, T ~ T c , for £=0.065 and also, 
albeit more moderately, for £=0.10. At these two compositions no Ts is observed and the sys tem remain tetragonal 
and paramagnetic down to T=0 K [89, '93J194]. Data in the SC state will be discussed in Sec. VII 

The Raman conductivity being independent of temperature for uj >600 cm -1 in the tetragonal phase, the quantity 
^dynamic can extracted in both symmetries by first extrapolating the data from the lowest energy measured 
(9 cm -1 ~ 1 meV) down to 0 cm -1 , and then integrating up to 600 cm 1 . The symmetry dependence of this 
quantity is shown in Fig. K3[c) for the parent compound Bal22. While the values of ^y narnic are very similar at 
high temperature in B\ g ana B^ g symmetries, the enhancement of the susceptibility upon cooling is only seen in B\ g 
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FIG. 6. Raman responses x" and Raman conductivity x" / u in the B\ g (a) and B^g (b) symmetry of BaFe 2 As 2 - The sharp lines 
superimposed on the electronic continuum are due to Raman active optical phonon excitations (c) Temperature dependence of 
the extracted static nematic susceptibilities in the two symmetry channels using Kramers Kronig relation (see text). 


symmetry clearly indicating an instability towards a charge nematic order with B\ g or x 2 -y 2 symmetry. 

The Co doping and temperature dependencies of ^^™ arnic are summarized in the color plot shown in Fig. J^b). 
The plot shows the clear correlation between the maximum of ^ d y™ arnic an d T s . It also highlights the persistence of 
significant nematic fluctuations over a wide range of Co doping covering most of the superconducting dome. 

As shown in Fig. [8ja) for each Co composition the nematic susceptibility can be well fitted using a Curie-Weiss 
law (see Eq. 16 in Sec. IIB and Eq. 42 in Sec. Ill B): 


Xb 


dynamic _ 


1 9 


= B + 


C 


T-T 0 


(55) 


where C = is a temperature independent constant. 

Here B is a constant which describes the temperature and symmetry independent part of the susceptibility, i.e. 
the non-critical part of the Raman response, and To is the Curie-Weiss temperature corresponding to the lattice- 
free electronic nematic transition temperature as defined in Sec. EH The extracted values of To follow qualitatively 
the doping dependence of Ts and extrapolate to zero slightly below optimal doping, at the critical doping x—x c ~ 
0.055 (Fig. J8^b)). As displayed in Fig. [8] for x=0.02 the temperature dependence of the nematic susceptibility does 
not depend appreciably on the incident photon energy ujl used for the Raman experiment. The insensitivity to ujl 
suggests that the resonant terms in the Raman vertex, if present, do not alter appreciably the temperature behavior 
of the Big response above Ts, justifying the use of the effective mass approximation as done in Eq. 


Eq. [52] of HID 


33 of III A and in 


It is important to note that the To values are systematically at least 40 K below the actual thermodynamic 
structural transition temperature Ts- This difference is a natural consequence of the absence of contribution of the 


lattice to the extracted susceptibility as discussed in |III C| To represents the bare nematic transition temperature of 
the purely electronic system while in the presence of the lattice, the actual transition temperature Ts is moved to 
higher temperature due to the finite coupling Aq between the electronic and lattice sub-systems (see Eq. 47). In 


that case the observed divergence of the nematic susceptibility is therefore cut-off by the structural transition which 
occurs at Ts > To. The difference between the two temperatures is a measure of the charge-lattice coupling energy 
x l/( c s r ~o9o)- This quantity, which is accessible both from Raman and from elastic constant measurements, is shown 
as function of Co doping in Fig. [8jd). 
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FIG. 7. (a) Temperature dependence of the Raman conductivity x" / u in the B\ g symmetry of Ba(Fei_ ;r Co ;r As )2 (Co-Bal22) 
for different Co doping concentrations x. (b) Color plot of the nematic susceptibility ^/™ armc as function of temperature and 
Co doping. 
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FIG. 8. (a) Curie-Weiss fits of the nematic susceptibility in the tetragonal phase of Co-Bal22. Structural transition temper¬ 
atures Ts are marked by arrows for x < 0.045. (b) Phase diagram showing the Co doping dependence of the bare electronic 
nematic Curie-Weiss temperature To (green squares), (c) Temperature dependence of the nematic susceptibility deduced from 
measurements at two different incident photon energies hr x=0.02 (T,s=110 K). (d) Co doping dependence of the difference be¬ 
tween the measured structural transition temperature T s and To. The open squares correspond to the same quantity extracted 
from a Curie-Weiss analysis of the shear modulus C s by Yoshizawa et al. [53] 
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B. Comparison with elastoresistivity and elastic measurements 


The observed Curie-Weiss like enhancement of the charge nematic susceptibility is qualitatively consistent with 
both elastoresistivity and elastic modulus measurements performed on Co-Bal22 [52H541 62 . The presence of orbital 
fluctuations above Ts was also inferred from point-contact spectroscopy measurements in the same system [95] • In 
the case of elastoresistivity measurements the purely electronic nematic susceptibility could be obtained by strain 
dependent measurements of the resistivity anisotropy in the tetragonal phase [62]. The extracted divergence was 
taken as evidence for a electronic driven structural transition. Since the nematic susceptibility extracted from Raman 
measurements is free from lattice effects, they also confirm the present of diverging electronic nematic degrees of 
freedom in the tetragonal phase of undoped and electron doped Bal22. We note however that the Curie-Weiss 
temperatures extracted from transport measurements are significantly higher than the ones obtained from Raman 
measurements on samples with similar Co doping. The discrepancy could be due to different couplings to electronic 
nematic degrees of freedom. Indeed, while Raman scattering couples to charge nematic degrees of freedom, the nematic 
component of the elastoresistivity tensor is a more complex quantity. On the one hand it could be a measure of Drude 
weight anisotropy arising from the sensitivity of the electronic structure to generate anisotropy in the presence of 
external strain [961. On the other hand, elastoresistivity is also sensitive to anisotropy in transport lifetimes of the 
carriers that can arise from scattering with the spin fluctuations or the impurities | 371 l 38 llMll 97 mnn]. Deviations 
from Curie-Weiss behavior have also been recently reported in several Fe SC near optimal doping by elastoresistivity 
measurements m ■ They have been interpreted as due to random field disorder that is a relevant perturbation near 
a quantum critical nematic transition. 

Elastic modulus measurements (discussed in a separate contribution in this issue [55]) can also be fitted with a Curie- 
Weiss like temperature dependencies. While Young’s modulus measurements from three point bending technique are 
now available for different Fe SC systems and dopings [54, 102 ] . direct measurements of the shear modulus Cs from 
ultrasound velocity measurements are only available for Co-Bal22 [21! [32 [53]. In both types of measurements a Curie- 
Weiss like softening is observed in Co-Bal22. The extracted Curie-Weiss temperatures are very close to Ts, which is 
consistent with a second order structural phase transition. We note that recent neutron scattering measurements in 
parent 122 compounds also observe a clear softening of the associated transverse acoustic phonon at low wavevectors 
[ill. Assuming a simple Landau type approach of the coupling between shear modulus and electronic nematic 


degrees of freedom identical to the one presented in |III C[ the bare electronic nematic transition temperature could 

-' n T—T cw cwi 

also be extracted from these elastic measurements using the relation Cs = Cg( t-t q )? where Tg ~ Ts and To is 
the purely electronic nematic transition temperature. The To values obtained from shear modulus measurements in 
Co-Bal22 agree remarkably well with the ones extracted from Raman measurements 1531E5I m- For both Raman 
and shear modulus measurements Ts~Tq increases mildly with Co doping (see Fig. |8jd)), indicating a possible incre ase 
of the electron-lattice coupling energy scale c o^ ogo (see Eq. 47) upon electron doping. We will come back in VIC for 
a more quantitative comparison between elastic and Raman measurements. 


C. Impact of disorder: Sr(Fei_ £C Co x As)2 

As already stated above, mechanisms of anisotropic scattering are a possible source of the transport anisotropies [97b 
HQO], associated with lifetime effects, that have been reported in both the orthorhombic and the tetragonal phases 
under applied strains or stress. Indeed transport measurements performed on annealed crystals appear to show a 
much reduced anisotropy in the orthorhombic phase, hinting at a key role of scattering mechanisms in the observed 
nematicity [38] • STM measurements have further shown the nucleation of nematic nano-domains around defect sites 
in the orthorhombic, and possibly even in the tetragonal phase [To5irm6| . These anisotropic impurity states were 
postulated to be responsible for the observed transport anisotropies because they are expected to act as strongly 
anisotropic scatterers eg mm ms]. Recent NMR measurements also suggest the presence of short range but static 
nematic order above Ts . This short range order could be due to pinning by impurities or to the presence micro-strains 
m and could possibly explain the onset of anisotropy observed at T > Ts in magnetic torque measurements mni. 

The above measurements have raised the question of the intrinsic nature of nematicity in Fe SC. A direct comparison 
between the Raman measurements in Bal22 (Ts= 138K) and Co-Srl22 (x=0.04, T£=137K) allows a direct assessment 
of the possible role of disorder in the emergence of nematic fluctuations above Ts m- Indeed despite having 
essentially identical Ts, resistivity measurements shown in Fig. [DJa) indicate a much higher residual resistivity ratio 
(RRR) in Co-Srl22 likely caused by the insertion of Co in the FeAs plane. Despite an order of magnitude difference 
in RRR, the extracted nematic susceptibilities shows extremely similar temperature dependencies in both systems 
(see Fig[9](b)). In particular the extracted Curie-Weiss temperatures Tq agree within ±5 K, clearly demonstrating a 
relative insensibility of the nematic fluctuations and their associated diverging susceptibility to disorder. The Raman 
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(a) 



(b) 



FIG. 9. (a) Temperature dependence of the resistivity of Bal22 (red, RRR~ 9) and Co-Srl22 (x=0.04, blue, RRR^l). Unpub¬ 
lished data courtesy of F. Rullier-Albenque. (b) Corresponding nematic susceptibilities extracted from Raman measurement 
in the B\ g symmetry uni- The lines are Curie-Weiss fits of the data points. 


measurements are consistent with recent elasto-resistivity measurements which show that while the magnitude of the 
observed transport anisotropies might be disorder dependent, the diverging behavior of the extracted susceptibility in 
the tetragonal state is not m • This conclusion is supported by a recent strain-dependent optical conductivity study 
which suggests that the strain induced transport anisotropy observed above T$ is not due to anisotropic scattering 
rate, but rather to anisotropic Drude weight E2- It is likely that disorder helps revealing an underlying intrinsic 
nematicity in transport and local probe measurements. The magnitude of the observed anisotropies is thus not 
necessarily the right quantity to assess the intrinsic nematicity of a given Fe SC system. 


VI. DYNAMICAL ASPECTS OF THE CHARGE NEMATIC RESPONSE 

A. Quasi-elastic peak 

We now go beyond the static properties of the nematic susceptibility discussed in the previous section, and analyze 
the frequency dependence of the nematic fluctuation spectrum revealed by the Raman measurements in B\ g symmetry. 
The temperature and symmetry dependencies of the Raman response in Bal22 (Fig. |6| suggest the presence of two 
contributions. The first one, broad in energy and weakly temperature dependent, is seen in both symmetries and 
dominates the Raman responses at high temperature. It can be assigned to intra and/or interband quasiparticles 
excitations which do not become critical upon cooling. This contribution shows only a weak suppression in the 
orthorhombic phase, which is linked with the simultaneous opening of the SDW gap. The second contribution, only 
present in B\ g symmetry, is strongly temperature dependent and is responsible for the strong enhancement of the 
static nematic susceptibility as discussed above. This contribution, obtained by subtracting the first contribution, is 
shown as a function of Co doping in Fig. [To|^a). The observed dynamical nematic fluctuations are quasi-elastic in 
nature and can be well reproduced by a damped Lorentzian lineshape of width Y and amplitude A in agreement with 
the form obtained in the theory part of this review (Eq. |40| ): 


Xqep 


= A- 


ujY 




r 2 ' 


( 56 ) 


The temperature dependence of the width of the quasi-elastic peak (QEP) T is shown in Fig. [lQ|^b). We find that 
T softens considerably upon approaching T$ for x < 0.045. This is because Y is approximately the single particle 
lifetime renormalized by (a/£ n ) 2 , where a is the unit cell length and £ n is the nematic correlation length, and the 
latter starts to increase as the system approaches the nematic instability with lowering temperature. Assuming a 
Curie-Weiss T dependence of the nematic susceptibility, and sufficiently close to To, the temperature dependence of 
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FIG. 10. (a) Evolution of the quasi-elastic peak contribution (QEP) in the B\ g symmetry as a function of temperature and 
Co doping. The contribution was extracted by subtracting the raw Bi g response by the one in B^ g symmetry for each doping. 
Fits using a damped Lorentzian are shown for each spectra, (b) Temperature evolution of width T of the quasi-elastic peak, 
(c) Temperature dependence of the inverse of the area A -1 of the QEP. (d) Co doping evolution of the mean-field nematic 
transition To extracted by three different methods. 


T is expected to vanish linearly with a zero intercept at To, the bare electronic nematic transition temperature (see 


Eq. 41 in III B ): 


T = 




( r o + c 3 )t 


(X (T-To). 


(57) 


For all Co compositions F is found to decrease linearly between T$ and up to at least 40K above T^. As in the case 
of the static susceptibility, the softening of F is not complete because of lattice effects which move the transition 
temperature to a higher temperature T^. 

The above analysis also allows us to extract the temperature dependence of the QEP area A , which is directly 
proportional to the diverging part of the nematic susceptibility i.e. without the background non-singular contribution 
B in Eq. |55| 


dynamic,QEP 
X Bl g 


[ Xqep( u ) ~ A 

Jo W 


(58) 


10 


c), the temperature dependence of A 1 is fairly linear, as expected for a Curie-Weiss behavior, 


As shown in Fig. 

with A -1 ^ (T — To). Extrapolation of the measured temperature dependencies of F and A~ 


_1 in the tetragonal phase 

down to zero temperature provide a measure of the scale To, which is in good agreement with the values obtained 
from the Curie-Weiss fits of the static susceptibility described in the previous section (see Fig. p!0[d)). 
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FIG. 11. (a) Temperature dependence of the slope of the Raman response at low energy Si. The solid lines are fits using 
Si = a' + ( T _ c Tq ) 2 • (b) Temperature dependence of the product FA where T is the effective one-particle scattering rate and A 
the diverging part of the nematic susceptibility (see text) 


B. Adequacy of RPA theory of the nematic transition 


Both the frequency and the temperature dependencies of the nematic susceptibility are in broad agreement with 
the expectations from the simple RPA description of a charge nematic transition outlined in the theory section of 
this review. The fact that the diverging nematic susceptibility is interrupted by the structural transition is a direct 
consequence of the lack of coupling to the lattice of a nematic susceptibility extracted from Raman in the dynamic 
limit. Both the Curie-Weiss behavior of the susceptibility and the linear decrease with T of the effective one particle 
scattering rate T are also in agreement with a simple mean-field description of a charge nematic instability described 
in |III B| The RPA type picture can be further tested by looking at the low energy slope of the B\ g Raman response, 
Si, which is predicted to vary as (T — To) -2 (see Eq. 43). As illustrated in Fig. [ll[b) the quantity S\ follows nicely 
the expected behavior for all Co dopings, even when forcing Tq to the value previously obtained from the Curie-Weiss 
fits of the susceptibility. 

A further test of the RPA picture can be made by looking at the temperature dependence of the product F.A. 
Identifying Eq. [40]with Eq. [56] this product should become almost temperature independent upon approaching Tq: 
FA rsj ^ _ 4 _ This is indeed what is found over a wide range of Co composition as shown in Fig. 11 [b). The 


decrease in — from x = 0 to x=0.045 can come from different sources. It could be linked to the decrease of the bare 

9or 

quasiparticle scattering rate \ as suggested by transport measurements in Co-Bal 22 [ 12 , 89], but it could also be due 
to an increase of the coupling constant go upon Co-doping. Overall there is, thus, a remarkable consistency between 
the behavior of the nematic B\ g Raman response in Co-Bal 22 and the theoretical expectations from a simple RPA 
theory of the nematic transition. 


C. A consistent picture with shear modulus measurements 


The Raman measurements indicate the presence of genuine electronic nematic fluctuations in the tetragonal phase 
of Co-Bal 22 . These fluctuations are expected to soften the lattice and trigger the structural transition [23, 76] . The 
availability of both Raman and shear modulus data in Co-Bal 22 allows us to draw a more quantitative picture linking 
electronic and lattice degrees of freedom. Indeed by symmetry the charge nematic order parameter O n probed by 
Raman and the orthorhombic lattice distortion eo are linearly coupled (Eq. [46]) . Remembering that Raman probes 
the bare charge nematic susceptibility, and using Eqs. (35) and (48), we have a simple relationship between x ( ^™ arnic 
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(a) (b) 




FIG. 12. (a) Mean-field model of the linear coupling between the inverse lattice (or shear modulus) and electronic nematic 
susceptibilities. The bare shear modulus C%, initially independent of temperature, is softened by the linear coupling to 
a diverging nematic electronic susceptibility -^^ namzc causing a phase transition at Ts > T 0 (b) Comparison between the 
experimental temperature dependence of the shear modulus Cs [53] for Co-Bal22 (x=0 and x=0.038) and the theoretical 
temperature dependence using the charge nematic susceptibility deduced from Raman measurements (x=0 and x=0.045, see 
text for details). To account for the slight mismatch in the structural transition temperatures Ts of the crystal used in the two 
experiments, the temperature axis was set in units of Ts. 


and Cs : 


n ^0 dynamic 

Cs = C S~ 


(59) 


where Ao is the charge lattice coupling constant introduced in |IIIC| and t\ is the nearest-neighbor hopping parameter 
introduced after Eq. (33). Assuming a weakly temperature dependent C$ the renormalized shear modulus can 


become soft via the linear coupling to the charge nematic susceptibility Xb narnic causing a second order structural 
phase transition at Ts defined as Cs{T = Ts) = 0 (see Fig. 12ja)). Fig. |l 2 |fo) shows a comparison between the 
shear modulus data of Yoshizawa et al. [53] and the theoretical shear modulus expected from the Raman data using 
equation 59 using Aq as the only free parameter. The quantitative agreement between shear modulus data and the 


theoretical expectation for both undoped and Co doped Bal 22 demonstrates that the charge nematic fluctuations 
observed by Raman scattering can fully account for the structural softening observed in Co-Bal 22 . It also provides a 
further evidence that the structural distortion observed in Fe SC is electronic driven. We stress that such quantitative 
agreement cannot be reached if we assume that Raman scattering couples to the fully dressed nematic susceptibility 
i.e. including lattice effects. The realization that Raman probes the bare nematic susceptibility is therefore crucial in 
reaching a consistent picture. 

One question left open in our discussion until now is the role of spin degrees of freedom in the nematicity observed 
in Fe SC. Anisotropic spin fluctuations associated with the stripe AF instability at or below Ts of many Fe SC have 
been invoked early on as a main driver for nematicity [SB HQ. Naively one might conclude from the consistency 
between shear modulus and Raman measurements that only charge / orbital degrees of freedom as probed by Raman 
scattering are relevant in driving the structural transition. This would however be incorrect, since it is very plausible 
that the diverging charge nematic susceptibility is itself enhanced or even driven by spin-nematic fluctuations via 
som-charee couDline’ [32l 111344118] . In fact it has been areued that a similar consistencv as in Fief. &) can be 
reached by comparing NMR measurements of spin fluctuations at Qsdw=(^^) /(tt, 0) and shear modulus data [ 119] . 
Disentangling between the electronic degrees of freedom behind nematicity thus remains a formidable task in 122 
systems m- Other systems like FeSe with no magnetic transition [16] [18j [82] 11021112044122] , or hole doped Bal 22 
with a still mysterious C 4 magnetic phase (12344126] may prove fruitful playgrounds to answer this question. It is also 
possible that, like the structure of the SC gap, the respective weights of various electronic degrees of freedom behind 
nematicity are not universal and depend on the Fe SC systems considered. 
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D. Digression: Raman fingerprints of nematicity in cuprates 

The presence of nematic correlations as possible precursor of uni-axial modulated charge / spin ordered states is 
also intensively discussed in the context of cuprates ]127H133j . It is therefore interesting to discuss here earlier Raman 
data on Lai-^Sr^CuC^ (La-214) and YBa 2 Cu 306+<5 (Y-123) which are somewhat reminiscent of the results discussed 
here, albeit with some key differences [65| f73l [13414136] . In underdoped La-214 a Raman QEP was observed to grow 
upon cooling over a relatively broad doping range [65] . The symmetry of the QEP was found to switch from B 23 to 
Bi^ channels at the critical doping p=0.05, indicating a 45 degrees rotation of the nematic fluctuations: along the 
Cu-0 bond at moderate doping and at 45 degrees of the Cu-0 at low doping. This is in agreement with neutron 
scattering studies where stripe-like spin modulations were also found to rotate upon increasing doping M- In La- 
214 there appears thus to be an intimate link between nematic (g=0) fluctuations observed by Raman and finite q 
stripe-like charge / spin fluctuations [ 1341 , 138{ . This phenomenology is similar to the spin-nematic scenario for the 
structural transition in Fe SC discussed above. In contrast with underdoped Co-Bal 22 however, the growth of the 
QEP in La-214 and its softening saturate at low temperature and therefore no true static long range order sets in 
|65 . Besides the growth of the QEP in La-214 is accompanied by the opening of a pseudogap at intermediate energies 
upon lowering temperature whose relationship with nematicity is still unclear. 

In other cuprates like Y-123, Bi 2 Sr 2 CaCu 208+<5 and HgBa 2 Cu 04 +( 5 evidence for nematic fluctuations in the Raman 
spectra are more elusive as no clear QEP is observed for moderately underdoped systems in any symmetry fTTiRirnffli 
Eoj. Notably in moderately underdoped Y-123 where CDW fluctuations are clearly observed in X-ray measurements 
and where transport anisotropies have been reported, no QEP is observed in the nematic B\ g and B 23 symmetries 
]135i 1139] . It has been recently suggested that the opening of the pseudogap may suppress the nematic QEP in 
Big symmetry in moderately underdoped cuprates [136 . We also note that in Y-123, whether the incipient order is 
nematic or bi-axial is still controversial jl41H149] . By contrast a clear QEP emerges in B 2 g symmetry for strongly 
underdoped composition (p < 0.05) [135] . In this regime however neutron scattering data indicate the presence of a 
uni-axial magnetic correlations along the Cu-0 bond rather the Cu-0 diagonal as suggested by the B^g symmetry 
of the Raman QEP [150 : . Here again the connection between Raman results and other probe appears problematic. 
The role of nematicity in cuprates and its link to the pseudogap is therefore largely unsettled as many competing and 
fluctuating spin and/or charge ordered phases are in close proximity. 


VII. NEMATICITY AND SUPERCONDUCTIVITY 

As of date, the study of nematicity in the Fe SC has concentrated mostly on the role of the various electronic 
nematic degrees of freedom in driving the structural transition in the normal metallic state. Relatively less attention 
has been given to study the nature of these fluctuations, if at all they exist, in the superconducting phase. In the 
following we show that, in a fully-gapped superconductor, soft nematic fluctuations lead to a new collective mode 
whose signature is a finite frequency resonance peak in the electronic Raman response. This behavior is qualitatively 
different from the quasi-elastic peak that develops when the nematic transition is approached in the normal state. 

A second motivation to study nematic fluctuations in the superconducting phase comes from the possibility that 
such fluctuations might affect the pairing mechanism and the superconducting transition. Theoretically, it has been 
argued that nematic fluctuations can play a role in the interplay between s and d-wave pairing states which are 
likely nearly degenerate in several Fe SC systems [15111158] . Concerning the effect on pairing, recent theoretical 
works [T59lirai| suggest that nematic fluctuations enhance superconductivity irrespective of the gap symmetry, at 
least for weak coupling superconductors. Therefore, it is possible that nematic fluctuations provide a complementary 
contribution to the pairing interaction on top of the (7r,0) magnetic fluctuations that promote pairing symmetry. 

On the experimental side, both X-ray measurements of the orthorhombicity and elastic modulus measurements 
from ultrasound velocity indicate competition between nematic and superconducting order parameters, which is in 
broad agreement with the expectations from simple Ginzburg-Landau-type theoretical arguments [23l l52l fTHTl [162] , 
However, these studies did not address the question whether the nematic/orthorhombic transition stays second order 
(or weakly first order) in the superconducting phase [163^ . and if yes, how the properties of the accompanying nematic 
fluctuations differ from those in the normal phase. In fact, experimental evidence for the existence of nematic 
fluctuations in the superconducting phase was lacking until recently. Consequently, it is notable that a recent study 
has identified resonance peaks in the electronic Raman response of Co-doped Bal 22 and Nalll as evidence of the 
existence of nematic fluctuations in the superconducting phase mu- 

With the above motivations, in this section we briefly discuss the theory of nematic fluctuations in a SC phase, and 
we compare it with the experimental findings from Raman response measurements in the SC state of electron doped 
FeAs. We finish by comparing the properties of the nematic resonance mode with those of other collective modes that 
have been postulated to exist in the context of superconductivity in the Fe SC. 
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A. Theory: Nematic resonance near a nematic quantum critical point 


In the superconducting phase the nematic fluctuations can be modeled by the same RPA approach as the one used 
to describe the QEP in the normal state data. As in the normal state, we assume that the nematicity is driven 
by electronic interaction of the form given by Eq. ft with the interaction constant go replaced by g , its value 
in the superconducting phase. Within RPA the interacting nematic susceptibility probed in B\ g symmetry in the 
superconducting state can be expressed in term of the bare superconducting response and g. It is given by (see Eq. 
[l0| and Eq. 34) 


XB lg M 


i -sn Bl9 M‘ 


(60) 


Mathematically it is clear that, since in the superconducting phase the Bogoliubov excitations are gapped, the structure 
of U B is quite different from that of a normal metal. For a fully gapped clean superconductor there are no fermionic 
excitations below 2A. Consequently, the imaginary part of the bare nematic response II^ is zero below 2A, and it 
is dominated by a pair-breaking peak at 2A. It is simple to check that, by Kramers-Kronig relation, this implies that 
the real part II^ diverges logarithmically upon approaching 2A from below. Therefore for positive nematic coupling 
g the RPA response develops a resonance below 2A at an energy Q r defined as 


i - g n' B jn r ) = o. (6i) 

This condition will always be met for a fully gapped superconductor in the (^-symmetric phase near the nematic 
instability. Thus, the presence of a nematic resonance in the Raman spectra is a generic property of the supercon¬ 
ducting excitation spectrum near a nematic quantum critical point. Physically, the opening of the gap shifts spectral 
weight from low frequencies uj ^ Y in the normal phase, where Y is the renormalized single-particle scattering rate 
defined in Eq. [4l] to higher frequencies uj ^ A in the superconducting phase. This shift of spectral weight transforms 
the quasi-elastic peak of the normal state into a resonance peak in the superconducting state. Close to the nematic 
quantum critical point the spectral weight of this nematic resonance quickly overwhelms that of the pair-breaking 
peak and grows as the critical point is approached as a function of the tuning parameter such as doping [164]. The 
presence of a full gap is however crucial for the existence of the resonance, as nodes will effectively wash out the 
divergence of the real part of the bare susceptibility II. For more details we invite the reader to consult Ref. m- 


B. Experiments: Fingerprints of a nematic resonance in superconducting Ba(Fei_ cc Co cc )2As2 


Raman scattering data provide a direct insight into this problem by following the evolution of the quasi-elastic peak 
upon crossing T c . Of particular interest are sample compositions which shows significant nematic fluctuations down 
to T c while staying in the tetragonal phase. This is illustrated for optimally doped Co-Bal22 (x=0.065) in fig. p~3[a). 
At this composition the low energy B\ g QEP observed in the normal state (see fig. 10) is strongly renormalized below 
T c with a suppression at low energy and the emergence of a well-defined peak at finite energy in the SC state. 

The SC peak is only observed in B\ g symmetry. It was initially interpreted as a Cooper pair-breaking peak 
associated to the creation of pairs of Bogoliubov quasiparticles in a BCS state [16511169] . Its energy is indeed close 
to twice the energy gap 2A detected by ARPES experiments at the electron pockets in similarly doped crystals 
[T70l [TO] , Its Co doping dependence however displays a striking departure from simple BCS-like expectations (see 
fig. [l3|(b)) . Coming from the tetragonal overdoped side (x=0.1, T c =20 K) the SC peak spectral weight is strongly 
enhanced upon approaching the boundary between the tetragonal and orthorhombic phases which located between 
x=0.065 and x=0.06 in Co-Bal22. In the nematic / SDW phase the peak spectral weight decreases and a much 
weaker peak emerges at a smaller energy [165] . 

The strong enhancement observed on the tetragonal side (x >0.06) of the phase diagram is at odds with the small 
changes in T c when going from x=0.1 (T c =20 K) to x=0.065 (T c =24.7 K). It is also inconsistent with the fact that 
the peak energy actually softens towards x=0.065 png pig. On the other hand the doping dependence of the 
peak spectral weight tracks the behavior of the nematic response above T c : it is stronger for the Co compositions 
when the nematic susceptibility is the strongest at T=T C , i.e. close to the nematic instability (see Fig. [I3[c)). This 
observation indicates that, as expected from the simple RPA picture discussed above, the same interaction responsible 
for the enhancement of the nematic susceptibility is also enhancing the SC peak spectral weight producing a nematic 
resonance close to the nematic critical point. We note that a similar enhancement of a B\ g SC peak spectral weight 
has also been recently observed in electron doped Co-Nalll upon approaching the nematic instability, indicating 
that the effect is possibly generic to electron doped FeAs systems m- This implies that the observed peak in 
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FIG. 13. (a) Evolution the B\ g Raman conductivity x" across T c for Co-Bal22 (^=0.065) [165] . (b) Co doping evolution 
of the SC Big Raman response (in blue). The responses just above T c are shown in black, (c) Evolution of the integrated SC 
spectral weight of the Raman response x" (hi blue) [165] as a function of Co doping. The corresponding nematic susceptibility 
Xdynamic ( n ) and well below T c (SC) are also shown in red. 


the superconducting B\ g channel is due to the nematic resonance rather than simply non-interacting pair-breaking 
physics. 

While the nematic resonance is essentially a delta function for a clean system, finite lifetime effects due to e.g. 
disorder are expected to broaden it. The broadening can be especially significant if the resonance is not fully detached 
from the 2A pair-breaking peak. It is interesting to note the peak observed in the experiments on Co-Bal22 are 
significantly broader than in Co-Nall 1: at least 5 meV (full width at half maximum) for the former while it can 
be as sharp as 1 meV for the latter mmm- This indicates that a fully detached and well-defined resonance 
is only present in Co-Nalll, while the nematic resonance manifests itself only as a enhanced pair-breaking peak in 
Co-Bal22. This may reflect a stronger coupling g in the case of Co-Nalll, but could also be due to a larger disorder 
due to the higher Co concentration of optimally doped Co-Bal22. 
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In principle the spectral weight of the resonance is expected to diverge and its energy to go to zero at the nematic 
quantum critical point x=x c . This corresponds to the situation where the nematic coupling reaches the critical value 
g c such that the Stoner criterion is satisfied ro=0 (see Eq. 15). However because the nematic Raman response does 
not couple to the lattice, this does not happen: the actual thermodynamic QCP occurs before the bare nematic critical 
point x c : xqcp > x c or qqcp < g c - In Co-Bal 22 x c can be estimated by extrapolating the doping dependence of 
To: x c rsj 0.055 [75l I164j . xqcp on the other hand lies at a slightly higher doping: between x=0.06 and x=0.065 (see 
Fig. [ 8 jb)) . Thereforejust like the divergence of the nematic susceptibility in the normal state is preempted by the 
structural transition which occurs at a higher temperature Ts > To, the complete softening of the resonance energy 
is also preempted by the nematic order which intervenes at a higher doping. This also explains the relatively mild 
softening of the mode energy seen in the experiments in Co-Bal 22 . 


C. Relationship with other superconducting collective modes 

The nematic resonance observed in the Raman spectrum of electron doped Fe SC bears a strong analogy with 
neutron resonance observed in several families of unconventional SC. In both cases they are due to residual electronic 
interactions in the particle-hole channel, which within RPA type pictures, trigger a well-defined excitonic-like pole in 
the charge and spin SC responses respectively. The nature of the residual interactions differs: magnetic and centered 
at finite qaf for the neutron resonance while they are centered at zero wavevector for the nematic resonance. Both 
resonances are connected to soft excitations associated to nearby quantum critical points: antiferromagnetic in one 
case and nematic in the other. It so happens that both critical points seem to exist in Fe SC and therefore both 
neutron and nematic resonance are present in their SC state. There is another key difference between both resonances: 
while the neutron resonance requires a sign-changing gap, the nematic resonance only requires the presence of a full 
gap. Because the presence of nodes is suspected in several Fe SC systems, it is likely that the nematic resonance will 
not be observed in several systems. The presence or not of a nematic resonance near the nematic instability is thus a 
powerful test for the presence of nodes for several systems including e.g. P-doped Bal 22 . 

Finally we should note that other collective modes have been predicted to occur in the SC state of Fe SC [17414178] . 
Among them the Bardasis Schrieffer (BS) mode |175l 111761 11791 Il80| , an in-gap electron bound state in a sub-leading 
pairing channel, is the most relevant for our discussion because of the near degeneracy of 5 and d-wave paring state 
in many Fe SC [4l 115211157] . Fingerprints of its presence have been reported in the Raman spectra of optimally hole 
doped K-Bal22 where several sharp modes are detected below T c |18ll 1182] . In principle both the nematic resonance 
and the BS mode can be simultaneously present in the SC Raman spectrum and may even couple [18311184] . However 
it is likely that their respective visibility will be strongly system dependent. While anisotropic gaps will suppress 
both collective modes efficiently, their presence is also linked to different criteria: the nematic resonance relies on the 
proximity of a nematic quantum critical point while the BS mode will occur when s and d-wave pairing states are 
nearly degenerate [TTHlUMj . In this respect we note that, from the point of view of nematicity, shear modulus and 
Raman scattering data shows much weaker nematic fluctuations near optimally hole doped K-Bal22 where a region of 
C 4 magnetic phase has also been observed [541 1186] . We therefore suspect that, by contrast to the electron doped side, 
the nematic resonance may not be present in hole-doped 122 systems at least near optimal doping. This illustrates 
the remarkable richness of Fe SC where the nature of the SC ground state and its collective excitation spectrum can 
vary strongly from one system to the other. 


VIII. CONCLUSIONS 

The purpose of this review was two-fold. First, we have provided a general theoretical framework within which 
one can understand electronic Raman scattering measurements near a nematic instability. Since Raman scattering 
measurement in the appropriate symmetry couples directly to the nematic order parameter involving the charge 
degrees of freedom, we have argued that it is one of the few direct probes of nematic fluctuations in a metallic system. 
However, contrary to thermodynamic probes like elastic constant measurements, the electronic Raman response probes 
the frequency-momentum dependent nematic susceptibility in its dynamical limit. An important consequence of this 
is that it is essentially blind to the lattice strains and the acoustic phonons. This, in turn, implies that Raman 
response is an unique probe of pure electronic nematic correlations. Experimental data on electron doped 122, and 
more recently 111 , Fe SC systems are consistent with this theoretical picture, and therefore provide unambiguous 
evidence of an enhanced nematic susceptibility of electronic origin. The comparison with shear modulus data confirm 
the absence of coupling between the electronic Raman response and the lattice strains. It also shows that the nematic 
fluctuations observed in the charge sector by Raman scattering can account for the observed lattice softening in 
electron doped 122 , providing further evidence that the structural phase transition is electronic driven. 
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However, despite the encouraging consistency between different probes of nematicity in the tetragonal phase, there 
remains several unsolved questions. The first one is related to the issue whether the electronic nematicity is driven 
by spin, orbital or charge degrees of freedom. In this respect FeSe is an interesting system because of the absence 
of magnetic ordering. While preliminary Raman measurements indicate the presence of charge nematic fluctuations 
above Ts in FeSe (122M187j , the nature of the magnetic fluctuations in this system remains to be clarified [1021118811189] . 
The second question is how generic are nematic fluctuations in Fe SC and whether or not they play a role in the 
SC pairing mechanism. Recent transport measurements advocate for the presence of underlying nematic quantum 
critical point near optimal doping for several systems and dopings. This will have to be confirmed by other probes 
and we believe electronic Raman scattering is ideally suited to answer this question. Indeed we have shown that the 
presence of a nematic resonance in the SC Raman response is a natural consequence of the proximity to a nematic 
quantum critical point, providing a smoking-gun experiment for its existence. 

From a more general perspective, the discovery of nematicity in Fe SC, first revealed via transport measurements, 
has stimulated intense efforts to design and understand experiments capable of probing nematic degrees of freedom. 
These efforts can now be capitalized to search for other systems where nematic degrees of freedom have been either 
predicted or indirectly observed. We have already mentioned cuprates and bi-layer ruthenates in the introduction, 
we can also add less correlated systems like bi-layer graphene [190H192] and even bismuth [Ml US]. The presence 
of nematic correlations can also be used to tune electronic orders via uniaxial stress or strain, just like magnetic 
field is routinely used to tune magnetic orders. Strain control of electronic phases has already been demonstrated 
in semiconductor based heterostructures like AlAs [ 1951 196]. In more correlated materials strain-effects have been 
demonstrated recently in Sr 2 Ru 04 where a strong T c enhancement under uni-axial strain was observed M- Up to 
now the microscopic origin of these kind of effects remain relatively poorly studied however. It is therefore likely that 
nematicity will remain a subject of intense research in the coming years. 
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